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SUMMARY 

V, 

V 

Recent  experimental  information  on  low-level  atmospheric  turbulence 
is  first  reviewed.  It  is  suggested  that  the  assumptions  of  homogeneity 
lii'.d  isotropy  customarily  adopted  for  high  altitudes  are  still  useful  in 
this  regime,  and  that  the  integral  scale  is  roughly  equal  to  9/10  of 
tie  altitude  up  to  about  1000  ft.  Next,  the  previously  published  theory 
of  the  ‘power-series  approximation’  as  applied  to  the  vertical  component 
of  t  ie  gust  is  extended  to  include  all  three  velocity  components  simul¬ 
taneously.  Fourteen  different  one-dimensional  input  power  spectra  and 
cross  spectra  are  found  of  which  only  five  are  important.  Of  these 
five,  only  one  is  a  cross-spectrum  involving  two  different  velocity 
components  (u  and  v)  .  Formulae  for  them  are  calculated  and  curves 
are  presented.  The  ‘gust  derivatives’  required  for  calculating  airplane 
response  are  defined  and  discussed,  and  the  most  important  ones  are 
shown  to  be  simply  the  negatives  of  classical  stability  derivatives. 
Methods  of  approach  for  calculating  the  remaining  ones  are  suggested. 
Finally,  it  is  shown  that  the  dispersion,  or  probable  error  of  position, 
is  fundamentally  different  when  the  controlled  variable  is  velocity  or 
heading  than  when  it  is  position. 
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riltOKV  OK  TIIK;  KLIGHT  ok  AIRPLANKS  in  isotropic  TCRBliLENCt  - 
KKVIKN  AND  EXTENSION 

B,  Etkin* 


1.  INTRODICTION  « 

The  ilielit  of  airplanes  ihroutth  turbulent  air  has  been  a  subiect  of  prime  concern 
to  aeronautical  engineers  since  the  beginning  of  flight  itself.  The  attendant  problems 
of  structural  integrity,  flying  qualities  and  performance  receive  continuing  study. 

The  application  of  statistical  methods  during  the  past  decade,  in  particular  the 
methods  of  power-spectral  analysis  and  the  theory  of  isotropic  turbulence,  have 
brought  about  a  significant  advance  in  our  understanding  of  these  problems. 

The  theoretical  approaches  to  analysis  fall  into  two  categories,  according  to  the 
manner  of  specifying  the  *unit’  element  of  the  gust.  The  first*  uses  a  ‘gust 

impulse'  as  the  basic  element,  as  shown  in  Figure  1.  References  1.  2  and  3  are 

representative  of  analyses  based  on  this  method.  The  second  approach uses  the 
elementary  spectral  (Fourier)  component  illustrated  in  Figure  2  as  the  basic  element. 
This  is  the  one  which  has  been  taken  in  References  4,  S  and  6.  It  should  be  emphasized 
that  there  is  no  fundamental  opposition  between  the  two  formulations;  both  can  in 

l9&d  to  'th?  £&!S6  results  uccur^c^  of  which  not  on  this  choice 

but  rather  on  the  details  of  the  approximations  subsequently  made  in  the  analysis. 

It  is  the  opinion  of  the  author  that  the  ..second  method  has  some  advantages,  viz: 

(1)  The  mathematical  formulation  is  simpler,  and  hence  easier  to  understand  and 
to  use; 

(2)  It  is  easier  to  separate  the  elements  of  the  theory  that  are  essentially  wing 
theory  from  those  that  are  essentially  the  representation  of  the  turbulence: 

(3)  By  using  the  power-series  approximation  of  References  5  and  6,  extended 
herein,  the  accumulated  knowledge  of  aerodynamics  embodied  in  stability  and 
flutter  derivatives  is  easily  incorporated; 

(4)  Approximations  involving  certain  parts  of  the  frequency  spectrum  are  easily 

Incorporated.  i 

This  report  presents  a  brief  review  of  the  information  on  atmospheric  turbulence 
in  Section  2.  It  follows  in  Section  3  with  a  semi-qualitative  description  of  the 
two  basic  methods  of  analysis  mentioned  above.  Section  4  contains  an  extension  and 
generalization  of  Reference  6  to  cover  the  case  of  simultaneous  inputs  of  all  three 
gust  components,  and  Section  5  presents  some  information  on  the  flight  path  of  a 
vehicle  flying -in  isotropic  turbulence. 


*  Professor  of  Aeronautical  Engineering,  University  of  Toronto,  Canada 

*  Due  to  H. W.  Llepmann,  Reference  1 
**  Due  to  H.  S.  Rlbner,  Reference  4 


2.  THE  STRUCTURE  OF  ATMOSPHERIC  TURBULENCE 


It  is  obvious  that  if  we  wish  to  study  flight  in  turbulent  air  theoret icr 1 ly  we 
must  know  enough  about  atmospheric  turbulence  to  construct  a  reasonable  mathematical 
model  of  it.  For  this  purpose,  the  atmosphere  close  to  the  ground  (in  the  boundary 
layer  produced  by  the  wind)  iieeds  to  be  consiuered  sepaiatcly  fi  uui  tiiat  h^giiei  up. 


2. 1  Outside  the  Boundary  Layer 


There  is  little  to  be  added  to  the  picture  of  turbulence  at  higher  altitudes  which 
has  alread)'  been  so  competently  given  by  Press  and  his  co-workers  at  the  N.A.S.A. 

(An  abbreviated  account  is  given  in  Reference  5).  In  short,  it  is  a  reasonable 
approximation  to  regard  high-level  turbulence  as  homogeneous  and  isotropic  in  patches 
of  limited  extent.  The  velocity  of  airplanes  is  normally  sufficiently  high  that  the 
tuibuient  field  within  one  of  these  patches  may  be  regarded  as  constant  during  the 
time  of  passage;  and  the  statistical  properties  of  the  input  to  the  airplane  are 
assumed  to  be  independent  of  the  response  of  the  airplane  itself,  i.e.  they  are  the 
same  as  would  be  obtained  in  rectilinear  translation  at  constant  speed.  (This  is 
not  to  say  that  the  response  of  the  airplane  is  neglected  in  calculating  the  aero¬ 
dynamic  forces.  The  forces  associated  with  motion  of  the  airplane  are  included  as 
usual.)  The  probability  distribution  of  the  intensity  a  of  the  turbulent  patches 
is  dependent  on  the  route,  season,  altitude,  etc.  The  one-dimensional  spectrum 
function  for  the  Idteial  cuuipuueiiL  of  Cue  turbulence  which  Is  now  widely  accepted  is 


e„(k,) 


cr^L  1  +  3k^ 

2rT  (1  +  k^)^ 


(1) 


According  to  the  theory  of  isotropic  turbulence’  the  above  is  derivable  from  the 
more  basic  energy  spectrum-function.  The  latter  is 

8  k 

E(k)  =  -cr^L  - 5  (2) 

■n  (k^  +  1)^ 


In  terms  of  E(k)  ,  the  one-dimensional  spectrum  is  calculated  from  the  relations 


and 


eij(ki) 


(kj,kj,k,)dkjdk3 


(3) 


E(k)  - 

= - 

477  k’  ”  ^  J 


where  is  the  Kronecker  delta.  When  Equation  (2)  is  substituted  into 

Equation  (4)  we  get 


(4) 


77 
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-k^ 


(k' 
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'I>jj  (k  j.kj.kj) 


(5) 
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and  Equation  (3)  becones.  for  the  particular  energ)  spectrun  adopted. 


2cr"L  -  kikj 

‘J  *  jj  (k^  T  1)^ 


dk  ^dk  j 


<6) 


Proa  Equation  (6)  we  nay  also  obtain  the  coopanion  to  Equation  (I),  i.e.  the 
longitudinal  one-diuensional  spectrum 


1 

■n  1  +  kj 


(7) 


The  cross'spectra  ,  0,,  .  6,,  are  all  zero,  since  for  these  cases  the  integrand 
of  Equation  (6)  is  antisyaaetrical  with  respect  to  one  or  both  of  k^  .  k^  . 

Unfortunately,  there  is  insufficient  information  available  on  the  scale  L  of  the 
turbulence  in  the  atmosphere.  The  value  L  =  1000  ft  has  been  assumed  by  Press  and 
others  to  be  reasonably  representative  but  much  more  experimental  information  is 
needed.  It  should  be  pointed  out  that  this  Is  a  very  important  parameter,  since  it 
miqr  exert  a  dominant  influence  on  the  energy  available  at  the  resonant  frequencies 
of  the  airplane.  This  effect  is  shown  in  Pigure  3,  taken  from  a  Douglas  Company  ‘ 
ieport*”.  Purthermore,  the  accuracy  of  the  power-series  method  (Sec. 5)  is  dependent 
on  the  ratio  of  airplane  size  to  turbulence  scale. 


2.2  Near  the  Ground 

At  low  levels,  the  turbulence  resembles  that  which  occurs  in  boundary  layers 
adjacent  to  rough  surfaces  and  is  strongly  affected  by  the  terrain.  The  scale  and 
intensity  both  vary  rapidly  with  height  above  the  ground,  and  in  general  the  field 
Is  neither  homogeneous  hor  isotropic.  A  number  of  measurements  have  recently  been 
reported  of  statistical  properties  of  low-level  turbulence*”"'^  from  which  two 
useful  general  conclusions  can  be  drawn.  The  first  is  that  Equations  (1)  and  (7) 
are  fair  approximations  to  the  lateral  and  longitudinal  one-dimensional  spectrum 
functions.  The  second  is  that  the  scale  factor  L  in  these  equations,  up  to 
1000  ft  altitude,  nay  be  approximated  roughly  by 

L  =  0.9  b  (8) 

where  b  is  the  altitude.  The  evidence  for  these  conclusions  is  given  in  Pigures  4 
to  6.  Figures  4  and  5  show  comparisons  made  in  the  USAF-supported  Douglas  study 
between  measured  spectra,  and  those  given  by  the  equations.  The  agreement  as  to 
shape  is  encouraging.  Pigure  6,  which  contains  more  detail  at  the  low  wave  numbers, 
is  another  comparison,  using  kj@(kj)  as  the  ordinate,  and  the  ratio  altitude/ 
wavelength  as  abscissa.  The  experimental  data  are'  those  of  Panofsky'*  and  the 
heavy  line  is  Equation  (1)  with  L  =  0.93  h  .  This  value  of  L  corresponds  to  a 
maximum  of  the  curve  at  b  -  0. 25  .  This  seems  to  give  the  shape  of  the  experimental 
curves  well  enough  at  heights  as  diverse  as  1  metre  and  300  metres.  No  Importance 
should  be  attached  to  the  actual  ordinates  of  the  curves  in  these  Figures,  since 
none  of  them  has  been  pornallzed,  and  there  are  wide  variations  in  cr  (which  is  the 
area  under  the  curve  when  plotted  to  linear  co-ordinate  scales);  only  the  shapes  are 
■  I,, 
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significant.  Panofsky  also  gives  a  seai-eapirieal  foraula  for  the  variation  of 
intensity  with  height  and  ground  roughness  under  unstable  aeteorological  conditions. 


This  is 

'j 

1 

a 

V 

-  0. 226 

1  rtfF  H 

- Q 

where 

V 

=  mean  wind  at  height  h 

“o 

=  characteristic  roughness  length. 

The  questicms  of  hosogeneity  and  isotropy  are  nore  troublesoae.  The  evidence  shows 
quite  clearly  that  low-level  turbulence  reflects  the  nature  of  the  terrain.  If  the 
latter  is  hoaogeneous  and  isotropic,  then  the  turbulence  will  be  closely  axisywaetric. 
i.e.  independent  of  rotation  about  a  vertical  axis,  and  honogeneous  with  respect  to 
translations  in  the  horizontal  plane.  However,  the  scale  and  intensity  in  general 
vary  with  height,  and  hence  the  turbulence  is  not  truly  isotropic  and  the  theory 
leading  to  the  tne-dinensional  spectnin  given  in  Equation  (6)  is  not  valid.  In 
spite  of  this,  there  would  seen  to  be  no  recourse,  in  the  present  state  of  the  subject, 
but  to  use  the  isotropic  nodel  for  the  low-level  case  as  well  as  for  high  altitudes. 
The  covlexity  of  the  problen  Is  even  t^en  quite  sufficient! 

-  ' .  ’-N 

Equation  (8)  indicates  that  we  nust  be  concerned  with  turbulence  having  scales  as 
snail  as  200-300  ft.  At  such  snail  scales,  the  variation  in  gust  velocity  over  the 
airplane  becones  important,  and  analytical  methods  of  sone  refinenent  and  couplexity 
are  indicated. 

3.  THE  T«0  BASIC  METHODS  OF  ANALYSIS 
3.1  The  'Inpulse'  Method 

Let  oXj  Xj  Xj  be  a  co-ordinate  system  so  chosen  that  the  mean  wind  in  it  is  zero, 
and  such  that  ox^  is  the  mean  flight  path.  Let  the  airplane  be  regarded  as  planar, 
so  that  only  the  distribution  of  atmospheric  motion  (u[  ,  u'  ,  up  over  the 
horizontal  plane  ox^  is  of  interest.  The  impulsive  gust  element  at  point 
(X,  ,  x^)  then  has  components 

uJdXjdx^  ,  u'di^dx^.  u'dXjdx^ 

of  which  we  consider  one  at  a  time  (as  for  example  in  Figure  1).  .  Now  let  the  air¬ 
plane  come  under  the  influence  of  the  gust  element  when  the  c.g.  is  at  position 
(4^, 0,0)  .  Then  a  typical  aerodynamic  force  or  moment  associated  with  it,  e.  g.  the  Z 
coag>onent  (the  negative  of  the  lift)  is 

upXjdX2h(Xj^  ^  ,  -  ^,Xj) 

where  h(Axj,X2)  is' the  response  function  for  a  unit-impulse  gust,  and  Is  zero 

for  Ax,  <  0  .  The  total  force  Z(x,  )  acting  on  the  airplane  Is  then  obtained 
'  c.g. 
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taOF  inteKratiiif  «ltb  respect  to  across  tbe  span  and  with  respect  to  froa 
•oD  to  CD  .  The  autocorrelation  of  Z<Xj^  ^  )  Is  next  obtained,  rix. 

2  \  ■ 

and  finallr  the  spectral  doisitjr  (ehich  is  the  quantitjr  soui^t)  is  obtained  bj  taking 
the  Fourier  Transfora  of  the  autocorrelation.  This  procedure  entails  soae  quite 
coaplicated  aatheaatics.  It  is  aorth  noting  that  the  basic  aerodynaalc  Infoiaation 
is  all  bound  up  in  the  function  h(Axj,x^)  . 

Thus  the  aethod  does  not  lend  itself  readilj  to  incorporating  aerodynaaic  inforaa- 
tion  (experiaental  or  theoretical)  ehich  is  in  the  fora  of  stability  or  flutter 
derivatives.  There  is  a  large  body  of  such  inforaation.  and  to  be  able  to  dran  on  it 
easily  Is  an  advantage.  Furtheraore,  ehen  ve  wish  to  extend  the  iwulse  aethod  to 
Include  the  three  velocity  covonents  slaultaneously.  the  coaplexity  is  further 
increased  by  the  preaoice  of  non-vanishing  teo-point  cross-correlations  between  the 
u[  .  u|  and  u'  coap(nents. 

3.2  The  TOarier  On^wneat*  Method  ^  „  » 

In  this  aethod  the  basic  eleaent  of  the  turbulent  velocity  field  is  a  wave  of 
shearing  notion,  described  by  the  expression 

e^°*dU(f^  (9) 

llie  corresptxidlng  distribution  of  downwasb  over  the  oXjX,  plane,  for  exaaple,  is 
shown  in  Figure  7.  Once  the  lift  and  other  relevant  sero^aaic  forces  or  aonents  . 
have  been  detemined  for  such  basic  velocity  fields,  the  foraalisn  for  writing  down 
the  spectra  of  the  Inputs  to  the  sirplsne  systea  is  quite  strsightforward.  However. 

In  Itself  this  step  does  not  aake  the  determination  of  the  basic  lift  eleaent  any 
easier.  It  replaces  the  problea  of  finding  biAx^.Xj)  with  that  of  finding  the 
periodic  lift  (or  other  force)  associated  with  a  running-wave  boundary  condition. 

In  fact,  the  latter  solutions  nay  be  constructed  hy  s  suitable  integration  of  tbe 
fomer.  Bxaaples  of  solutions  of  this  kind  of  wing-lift  problea  are  found  in 
References  14  and  15. 

3.2.1  The  Power-Series  Approximation 

A  slaplifying  approxiastlon  introduced  in  Reference  5  and  extended  in  Reference  6 
is  based  on  representing  the  gust-velocity  field  over  tbe  airplane  by  a  nodified 
Taylor  series.  It  was  shown  in  Reference  6  that  hy  keeping  tens  in  tbe  series  up 
to  the  second  order,  tbe  velocity  distribution  can  be  represented  adequately  for 
spectral  couponents  whose  wavelengths  on  the  two  axes  (X,  and  X^  .  Fig.  7)  exceed 
twice  tbe  corresponding  airplane  dimension  (length  or  span)*.  It  was  further  shown 
that  the  cut-off  frequency  obtained  by  excising  the  higher  wave  numbers  is  high  enough 
to  allow  inclusion  of  important  elastic  modes,  and  that  tbe  error  due  to  using  a 


*  It  Is  shown  later  (Sec. 4. 1.1)  that  the  wave-length  llaitation  is  actually  less  restrictive 
than  this. 
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truncated  spectrum  is  not  serious  provided  that  the  ratio  L/l  is  not  leas' than 
about  3.  The  value  of  1  for  a  largo  sveptaing  airplane  is  about  100  ft.  so  the 
turbulence  scale  L  aay  be  as  sbmII  as  300  ft  for  such  aircraft.  For  smaller 
machines.  L  nay  be  correspondingly  less. 

It  SuCuXu  tii&t  A  t  Buy  ficQututly  uoi  be  uecessary  iu  retain  the  second 

order  terms  in  the  power'series  development.  From  the  examples  shorn  in 
Reference  5  it  can  be  seen  that  cutting  off  the  spectrum  at  component  aave-lengths 
less  than  eight  times  the  wing  chord  and  wing  span  respectively  may  still  provide 
sufficiently  good  results  for  motion  in  the  rigid-body  modes.  This  requires  only 
that  the  zero  order  and  linear  terms  of  the  series  be  retained.  Furthermore,  it 
mill  be  seen  in  the  folloming  that  the  input  spectra  associated  mith  the  second- 
order  terms  are  very  small. 

./ 1  . 

4.  EXTENSION  •F\TiiE  POtEt-SERIES  APPMXiaATION 

In  References  S  and  $  only  th«  vertical  compeaent  of  the  turbulence  **  «') 
mas  considered  to  be  present.  Homever.  the  slMiltaneous  occurrence  of  all  three 
velocity  components  must  be  considered  for  a  complete  theory.  Thus  me  take  ms  the 
description  of  a  single  spectral  component  of  the  gust  field  the  Taylor  series 


u 


r 

1 


4“ 


If  \ 


i  =  1.2^3 
j.h  =  1.2 


(10) 


mhere  the  suamation  convention  for  repeated  suffices  is  implied.  The  subscript  o 
denotes  the  airplane  C.O..  l.e.  the  point  (Ugt.o.o)  .  Thus  n^  and  its  derivatives 
are  periodic,  mith  circular  frequency  n^u^  =  w  . 


As  in  Sectiew  3.1,  me  consider  the  airplane  to  be  a  planar  body,  so  that  only 
the  variations  in  the  plane  are  of  Interest  -  hence  the  restriction  of  J.  k  to 
1,  2  in  Equation  (10).  In  Reference  6  a  refinement  mas  Included  mblcb  improved  the 
fit  obtained  mith  this  approximation  to  the  actual  sinusoidal  velocity  distribution. 
The  refinement  mas  to  multiply  the  linear  terms  by  suitably  chosen  frequency- 
dependent  factors.  This  had  the  same  effect  on  the  Input  spectrum  functionB  ms 
mould  adding  certain  third  order  terms  to  Equation  (10).  Altbongdi  there  is  certainly 
some  gain  In  accuracy  achieved  thereby,  this  refinement  adds  undesirable  complexity, 
and  is  not  included  herein. 


The  point  of  vlem  taken  is  that  each  term  of  Equation  (10)  when  it  is  applied  to 
a  single  spectral  component)  represent^  a  periodic  relative  velocity  field  of  slvle 
form.  lAlch  results  in  periodic  aerodyil'amlc  forces  and  momenta.  These  are  expressed 
quite  generally  by  a  set  of  ‘pint  derlrstlves'  or  'gust  transfer  functions’  which 
are  analogous  to  (soaie  are  identical  tc)  the  faasillar  stability  and  flutter  deriva¬ 
tives  which  have  been  in  use  for  so  long. 


Consideration  of  the  syimietry  of  the  velocity  distributions  represented  by  the 
individual  terns  of  Equation  (lO)  permits  separation  of  the  associated  aerodynamic 
forces  and  moments  into  the  usual  longitudinal  and  lateral  groups.  The  following 


7 


utrix  equations  serve  to  define  the  'gust  derivatives*  (note  that  the  ;ust  velocities 
are  noe  denoted  by  u'  ,  v'  .  w')  ; 

{Pj}  =  [Ajl  {g,}  (11) 


{Pj>  =  [Ay]  igj} 


(12) 


where 
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(16) 
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where  etc.  It  will  be  seen  subsequent lir  that  the  input  spectra  associated 

with  m‘  are  negligible,  and  hence  that  tern  way  be  dropped. 
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In  the  above  expressions  {Fj}  and  {p^}  are  the  column  matrices  of  the  Laplace 
transforms  of  the  longitudinal  and  lateral  aerodynamic  forces  respectively,  {g^} 
and  {gj}  9.T6  the  matrices  of  the  Laplace  transforms  of  the  gust-velocity  inputs  for 
the  longitudinal  and  lateral  equations,  and  [A^]  and  [A,]  are  the  matrices  of 
'gust  transfer  functions’  defined  by  Equations  (11)  and  (12).  These  transfer  functions 
might  frequently  be  approximated  by  simple  derivatives,  e.g.  =  dY/dv"  (see 
Ref.  5,  Sec.  4. 16).  The  matrices  [a^]  and  [A^]  are  written  out  above  with  maximum 
generality,  in  which  case  there  are  a  total  of  80  transfer  functional  The  dashed 
lines  in  Equations  (15)  to  (18)  indicate  those  portions  of  the  matrices  (to  the  right 
of  the  line)  ahich  would  be  neglected  in  a  first-order  theory.  The  number  of  transfer 
functions  is  then  reduced  to  40.  If  only"control-fixed  conditions  are  of  interest,  a 
further  reduction  to  27  is  effected  by  dropping  all  the  H  terms.  Additional 
simplifications  of  the  sort  common  in  stability  and  control  work  might  frequently  be 
in  order;  for  example,  neglecting  the  X-force  equation  altogether  in  the  longitudinal 
equations  of  motion,  and  dropping  certain  transfer  functions  which  experience  or 
analysis  indicate  are  small. 


4. 1  Hie  One-Dimensional  Input  Spectra 

Since  the  'Inputs'  {g^}  and  {g^}  contain  more  than  one  element,  the  airplane 
system  is.  subjected  to  a  set  of  simultaneous  random  Inputs.  Figure  8  illustrates  the 
general  case,  with  inputs  x^(t)  ,  1  =  l  to  n  ,  output  y^Ct)  ,  and  transfer 

functions  G„i(s)  .  The  output  is  given  by 
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where 
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It  follon  froB  Equations  (11),  (12).  (21)  and  (24)  that  [O,^]  =  and 

[O,]  =  .  where  [q,1  and  [Oj]  are  the  aatrlces  of  the  oreraU  transfer 

faactl«i5.  for  the  teo  aeta  of  equations.  In  the  ateve  equations,  the 

quantities  Z,  .  etc. .  are  to  be  interpreted  as  transfer  functions,  1. e. 

dZ  Bz 

Z,^  =  —8  +  ~  .  etc.  (25) 

dw  v9 

The  equations  do  not  include  any  autoaatic  control  elesMnts,  but  the  addition  of  these 
'  in  particular  cases  Is  usually  fairly  straightfomard. 

.V  ’  /» 

■e  Hust  now  consider  the  input  spectrun  functions  which  occur  In  Equation  (19). 
These  are  the  cross-spectra  of  all  the  inputs  that  occur  in  (g^)  or  {g^}  ,  that  is, 
anong  the  velocity  cosponents  and  their  first  and  second  derivatives.  Mav  of  these 
cross-spectra  are  aero  by  virtue  of  the  fact  that  the  two  quantities  Involved  are 
uneorrelated  (see  after  Ekpi.  (7)).  However,  a  nunber  of  then  resiain,  and  these  aust  be 
calculated.  Let  the  spectrun  function  corresponding  to  any  pair  of  entries  in  {g^) 
or  {g,}  be  identified  fay  a  eorrespondinc  pair  of  subscripts.  For  exsople,  Bux^xy 
is  the  cross-spectrun  of  u'  and  v^  which  occur  In  {g^}  .  The  expression  for 
the  three-dinenslonal  cross-spectnia  of  two  scalar  coaponents  of  a  vector  of  the 
fora  given  by  Equation  (9)  is  given  by  Batchelor  (Ref.  9.  Eqn.2.5.  5).,  as 


.  X  dUj(i1)dD,(f5> 

#,.(10  =  lla  "  *  '  J  — 


The  cross-spectra  of  eleaents  containing  derivatives  can  be  written  down  directly 
froa  Equation  (26).  For  exaaple,  the  spectral  coaponent  of  u,  ,  froa  Equation  (9) 
is  given  by  the  x,  derivative  of  the  u,  coaponent,  viz. 


‘(in^dU,) 


whence  for  exaoiple 


I, 


Ibe  guisral  rule  Is  sssn  to  be  thst  for  each  derivative  with  respeot  to  the 
speotruB  function  is  aultlplied  by  .  Hie  plus  sign  is  jfor  derivatives 

of  the  second  subscript  velocity  cosponent  (v,^  in  Bqn.  (2B)).  um  the  sinus  sign 
is  for  derivatives  of  the  first  (u,)  .  The  difference  in  sign  ooburs  because  the 
eonjfugate  of  the  first  aaplitude  elesent  la  used  in  Equation  (26).  The  correspond* 
ing  one-disensional  spectnis  function,  continuing  with^the  same  exaanle.  is  then 


(cf.  Bqn. (3)) 
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^k3«33(k,.k3.k3>dk3dk3 

(b) 

W 

In  the  theory  presented  herein,  we  exclude  that  portion  of  the  spectnis  for  which 
>  D'  and  >  QJ  where  and  ()[  correspond  to  the  wave-length  Units  for 
idilch  the  power-series  approxlsatlon  is  valid.  It  sust  be  noted  that  some  of  the 
integrals  of  the  type  coutained  in  Equatlou  (2P)  are  divergent  when  the  Units  are 
infinite  andv'the  truncation  is  therefore  -  ssential.  and  not  a  natter  of  choice. 

The  expression  for  the  truncated  sps'^trun  is 

-kj  •« 

With  the  value  of  given  in  Equation  (5)  this  integral,  and  the  others  like  it 
which  occur  in  the  equations,  can  all  be  evaluated  quite  sisply.  The  integrand  in  the 
ssjprlty  of  eases  is  an  odd  function  of  one  Oi  both  of  kj  and  k,  and  for  these 
the  integral  is  zero.  Of  those  irtiicb  resaln,  sone  can  be  discarded  on  the  basis  of 
the  following  order-of-sagnltude  analysis. 


Ibe  general  fom  for  9(kj}  (apart  fros  sign)  is 


9(ki) 


(ikj)“(lk,)^»ljdk3 


(31) 


where  n  =  a  +  /3  ^  and  _  a  and  ^  are  the  orders  of  the  two  velocity  derivatives 
involved.  When  the  expression  for  given  by  Equation  (5)  is  inserted,  we  get 


e(ki) 


(k*  +  D* 


dk. 


(32) 
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Depending  on  the  values  of  i  and  j  .  the  integration  with  respect  to  k,  leads 
to  aero,  or  one  or  both  of  the  following  terns 

^  (kj +  k* +  !)-»/*  ;  ^  (kj  +  k*  +  n”/*  (33) 

Since  we  are  interested  in  values  of  kj  and  k,  up  to  about  100,  we  see  that  the 
■agnltude  of  6  is  characterized  br  the  nuabers 

L^"“(100)""'‘  or  l‘"“(100)“** 

of  which  the  larger  one  is  the  second. 

Thus  the  relative  values  of  0  with  ascending  n  are  character Ized  as  shown  in  the 
following  table: 

Tiable  I 


n 

0 

1 

2 

3 

Relative  6 

1 

0.1 

0.01 

0.001 

for  L  =  1000 

1 

0.5 

0.25 

0.125 

for  L  =  200 

On  the  strength  of  these  values,  and  noting  that  L  -  200  is  a  rather  snail  scale, 
we  may  neglect  all  cross-spectra  for  which  n  >  2  .  The  reaalning  non-zero  spectra 
(25  in  all)  have  been  calculated  and  are  given  below. 


Por  the  Longitudinal  Equations 


n  =  0 


n  =  1 


-uu 
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_  3^ [i  -  i 
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27rk'(l  +kj)[  \kV. 


0irw,  = 


3  a*  k,k^  r  k»  i  M /iiY  t 
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(34) 


(39) 


(36) 


(37) 
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(39) 
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(41) 

=  -  «UT„ 

(42) 
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(43) 
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(44) 

(45) 
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(46) 
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=  3<r*r£ji^/k'V  /kn  2/kq^^^k^.kn 
An\.\z  r  *  \\\k’ )  \kV  sVkV  k' -  k'J 

(47) 

For  Lateral  Eqpatione 

®irr 

^  3o-*L  1  k'f  kf  r  l/'k'\*l  ll 

2v  1  +k|  k'[l  +  kj  [*  slkOJ  *  s\ 

(48) 

0WyW,  =  -Qwyy  (58) 

The  spectra  given  above  are  plotted  against  for  several  values  of  k'  In 
Figures  9  to  22.  It  aay  be  noted  that  none  of  then  are  cov>lex  -  they  are  either  | 

real  or  pure  Inaginarles.  There  are  25  non-zero  poeer  spectra  and  cross-s^peetra  I 

I  listed  above.  IkniT  of  these  are  ettual  or  aerely  opposite  In  sign  to  others,  so  that  { 

I  there  are  cnly  fourteen  essentially  different  ones.  Of  these  fourteen,  three  are 

j  sero-order  (0uu  .  9rr  .  9ww}  four  are  first-order  •  ^vy  >  Bw,  .  6w«()  and 

j  the  renslnlng  seven  are  all  second-order.  Of  the  first-order  spectra,  only  one  Is  a 

I  cross-spectruB  Involving  two  different  velocity  covonents.  i.e.  0uvy  .  Hence  In  a 

I  first-order  theory,  this  reaains  as  the  only  cross-tern  between  velocity  co^Knents. 

I  and  if  it  is  neglected,  comlete  statistical  separation  of  the  response  to  the  three 

i  eonponents  of  the  turbulence  results. 

4.i.i  The  Wavelength  Limitation 

Exaaination  of  Figures  9  to  22  and  Table  1  shows  that  the  order  of  aagnltude  of 
the  spectrua  peak  is  given  by  .  Nov  if  the  basic  series  giving  the  velocity 

(Eqn.  10)  had  been  extended  to  Incldde  blgher-Ohdar  terns,  the  effect  would  siaply 

■  i 


IS 


have  bMn  to  add  additional  hlgiii«r>ordar  spactra  (a  >  3)  to  tbe  Hat  alraady 
calculatad.  It  la  avldent  that  those  hlgher-order  spectra  eoald  be  nexllglbli  sanll 
for  the  freqeenc;  range  k  <  1  and  for  the  acnle  L  >  100  .  In  the  range  1  <  k  <  100 
they  nould  nltlnatelr  hecone  large  as  n  Increnaed  Indefinitely.  Ihus  It  appenra 
that  the  spectra  presented  are  actually  valid  for  a  series  representation  of  the 
velocity  containing  terns  of  at  least  the  third,  and  profanhly  higher  order.  The 
eavelength  llnltatlons  nay  therefore  reasonahly  he  taken  as 


K>  =  b 


Aere  la  the  overall  length  of  the  nlag.  Hence 


\\ 


gj  =  2.-  =  3.- 


Fdr  exaaple.  If  L  =  1000  and  h  =  100  .  then  k^  =  20vr  =  62.8 


PinjHly.  it  nay  he  renarked  that  for  large  L  .  (l.e.  1000  ft)  the  seeond-order 
spectra  are  leas  at  aadlun  nave  nuabera  (k  1)  by  a  factor  of  order  10*  than  the 
aero  order  spectra.  Thus,  aniens  relatively  hUth  freoaeney  rw^onses  are  of  lattrsnt 
(e.g.  elastic  nodes)  the  aeoond-order  spectra  are  not  at  all  Isportaat. 


4.8  The  ^sst  Derivatives* 

Bguatlana  (17)  and  (18)  Indicate  that  the  general  second-order  theory,  ^mn  applied 
to  the  rigid-body  ■otico  of  an  airplane  nith  three  additional  control  degrees  of 
freedon.  involved  80  aerodynanlc  transfer  fanctlons  (nhleh  ne  have  tersed  ^gost 
derivatives').  Shoold  additional  elastic  degrees  of  freedon  of  the  airplane  he 
Included  (as  in  Ref. 6),  then  still  additional  gust  derivatives  eonld  be  regnired.  As 
has  already  been  nentloned.  boeever,  sobstantial  slvllfloations  can  be  nade  in  asny 
practical  analyses,  such  as  dropping  the  X-force  equation,  keeping  only  the  first- 
order  derivatives,  etc.  These  slnpllflcations  nost  aleays  be  deternlned  by  the 
particular  circuastances,  and  it  is  not  within  the  scope  of  this  paper  to  anticipate 
all  the  possibilities.  Neither  Is  It  within  its  scope  to  present  a  collnctloo  of  data 
OB  the  derivatives,  although  it  is  hoped  that  sane  research  at  the  Institute  of  Aero- 
physics  will  be  directed  to  that  end.  Nevertheless,  a  discnsslon  of  the  derivatives 
..is  giv«i  in  Sections  4.2.1  to  4.2.3. 


„  4.2.1  The  Zero-Order  Derivatives 

The  zero-order  derivatives,  which  are  the  aost  isportant  ones,  are  those  with 
respect  to  the  gust -velocity  conponents  theaselves,  e.g.  Il,»  .  .  etc.  They  are 

the  eleaents  of  the  first  two  coluans  of  [a,]  and  the  first  colnan  of  [A,]  . 


These  are  simply  the  aero-ivnaaic  transfer  functions  (staMlity  ie:  .vntives)  of 
classical  aerovtvnaaic  theory,  with  opposite  sign,  i.e. 

Z^>  =;  -Z,  .  etc.  (59) 

where  is  given  by  Equation  (25).  The  reversal  of  sign  is  because  w  is  the 

velocity  of  the  airplane  in  the  z-direction.  and  w'  is  the  velocity  of  the  air  in 
the  same  direction;  hence  the  relative  motion  is  given  by  (w  -  w')  .  This  group  of 
derivatives  embodies  the  major  aerodynamic  effects  of  gusty  air.  and  a  simplified 
calculation  in  which  all  others  are  neglected  would  still  be  of  considerable  value, 
especially  for  small  airplanes  in  large-scale  turbulence. 

‘i.?.’?  The  First-Order  Derivatives 

Colunins  3  to  5  of  [.Ajj  and  columns  2  to  4  of  [a^]  contain  the  elements  in  which 
there  aPi>ears  a  first-order  derivative  of  the  velocity  components.  These  describe 
the  influence  of  the  ‘gust  gradient*  on  the  airplane,  and  are  no  doubt  important  for 
large  airplanes,  especially  near  the  ground  in  small-scale  turbulence.  It  has  already 
been  shown*  that  the  derivatives  with  respect  to  w'  and  w^  '  are  identical  with  the 
classical  pitch  and  roll  stability  derivatives,  viz.: 


=  -M-  ,  etc. 

X 

and  (60) 


No  correspondingly  simple  interpretation  is  in  general  possible  for  the  velocity 
fields  associated  with  u.[  and  .  For  unswept  wings  of  high  aspect  ratio,  the 
derivative  u^  would  presumably  be  significant  only  in  introducing  a  relative  wind 
at  the  tail  different  from  that  at  the  c.g.,  i.e. 


Au 


rel 


t 


This  would  codify  the  tail  lift,  and  hence  the  lift  and  pitching  moment  of  the  air¬ 
plane.  as  expressed  in  the  derivatives  Zu^  and  Mu*  •  For  sweptback  high-aspect- 
ratio  wings  it  introduces  a  variable  (linear)  relative  wind  along  the  span,  which 
could  be  treated  by  a  suitably  modified  lifting-line  theory.  The  same  theoretical 
wing  problem  is  presented  by  the  velocity  field  associated  with  ,  with  the 
difference  that  the  spanwise  velocity  variation  associated  with  the  latter  exists  for 
all  wings,  whether  swept  or  not.  For  the  particular  case  of  a  straight  lifting  line, 
the  forces  corresponding  to  u'  are  just  those  given  by  the  classical  yaw-rate 
derivatives,  viz. : 


The  effects  of  the  linear  velocity  fields  associated  with  and  Vy  on  the 

contributions  of  the  vertical  tail  to  the  aerodynamic  forces  can  readily  be  estimated, 
since  they  merely  change  the  average  relative  sidewind  wind  at  the  tall  and  hence 
the  angle  of  attack  of  the  vertical  tail.  Their  effects  on  the  wings  are  rather  more 
involved,  v'  would  not  be  expected  to  be  of  much  Importance  for  unswept  wings,  but 
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for  swept  winis  sad  both  have  the  effect  of  aodifylns  the  wing  axi^le  of 
attack  distribution,  when  it  has  dihedral,  in  the  aanner  iil.xsti'a':  •’  in  Figure  23. 
Again,  for  hlgh«aspect-ratio  wings,  lifting-line  theorj  could  be  used  to  calculate 
these  effects  in  a  rather  straightforward  aanner.  For  score  general  cases,  lifting- 
surface  theory  would  have  to  be  eaployed. 

then  the  wing  is  swept  there  is,  in  addition  to  the  a  changes  described  above, 
the  iwportant  variation  of  the  aagnitude  of  the  coaponent  of  the  relative  wind  norml 
to  the  line  of  aerodynaaic  centres.  This  is  given  hy 

AVjj  =  +v'  sin  A  .  y  <  0 

and  the  distributions  of  AV^  associated  with  r'  and  have  exactly  the  snae 

fora  as  those  shown  for  Aa  in  Figure  23.  Thus  the  two  effects  will  be  additive  in 
producing  rolling  ocatent,  side  force  and  yawing  aoaent. 

4.2.3  I'he  Sec  Of  tJ -Order  Oerivat  iires 

Columns  5  to  10  of  [a^]  and  4  to  3  of  lAj]  contain  the  second-order  eleaents. 

By  virtue  of  the  a^uaption  cade  in  Section  4.1.  i.e.  neglecting  all  input  spectra 
having  n  >  2  .  one  coluan  of  these  derivatives  is  not  required.  That  is  the  seventh 
coluan  in  [a^]  ,  containing  derivatives  with  respect  to  w^  .  The  reason  for  this 
is  that  the  lowest -order  spectrua  function  which  contairus  the  input  w^  is  ^tw{j 
(n  =  2)  and  it  is  identically  zero.  Hence  this  particular  input  is  of  negligible 
iwportance  and  the  associated  derivatives  are  not  of  interest. 

Of  the  reaaining  derivatives,  those  invoiviag  w'^  and  w^  have  already  been 
dis5!i£isad  in  Reference  6  (using  a  different  nooenclature ) .  They  are  shown  to  give 
the  aerouynaaic  forces  resulting  frea  a  periodic  caabering  or  chordwise  bending  of 
the  wing  (Wjj)  and  a  flapping  or  spaawise  bending  •  Values  of  the  lift  and 

pitching  aoaent  on  a  two-diaensional  wing  in  incoapresslble  flow  are  given  there  for 

case.  The  calculation  of  forces  due  to  the  w^^  field  could  be  accoaplished 
hy  a  relatively  straightforward  application  of  the  appropriate  aethod  of  wing  theory. 

The  elevator  and  rudder  hlnge-ooaent  derivatives  contained  in  [A^]  and  [Aj] 
could  all  be  calculated  relatively  easily  on  the  assuaption  that  the  surface  in 
question  experiences  an  angle  of  attack  or  velocity  change  equal  to  that  at  the  nean 
aercdynaaic  centre  of  the  surface.  The  calculation  of  aileron  hinge-notsent  deriva¬ 
tives  (the  4th  row  of  [a^])  would  take  more  effort  except  when  slaple  strip  theory  i.s 
acceptable. 

Generally  speaking,  since  the  input  spectra  corresponding  to  n  =  2  are  relatively 
so  weak,  it  appears  that  rough  estimates  of  the  second-order  derivatives  will  serve 
well  enough  for  analysis.'  A  note  of  caution  must  be  sounded  in  this  connection, 
however,  when  elastic  modes  of  the  aircraft  are  Involved,  for  then  the  second-order 
terms  nay  be  more  important. 


Reproduced  From 
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5.  BISrUSItN  9€  THE  FLIMT  fATE 


nan  tb;  aircraft  is  flora  hy  a  huaan  or  autoaatic  pilot  so  as  to  travarse  a 
iqwclfiad  track  (e.g. ,  as  glTsn  hy  a  radio  beta),  at  a  specified  altitude  (e.g. ,  as^ 
giTe£  hy  as  altisctcr),  thts  the  ccstrollcd  Variables  sajr  be  cousiuercu  as 
(altitude)  aod  x,  (lateral  divlaceaeot).  These  will  be  randoa  rmrlablesi.  having 
■eea'SQoare  values  vfalch.  when  used  In  the  nonal  (Oaussian)  probability  distribution, 
give  the  probability  of  dlspersloB  of  the  aircraft  frea  the  desired  (rectilinear) 
flight  path.  In  a  hoaogeneous  isotropic  ataospbere  this  probability  functiOQ  applies 
equally  veil  to  all  portions  of  the  path.  Boeever,  afaen  the  controlled  variable  is 
a  velocity,  rather  than  e  displacenent,  the  sitimtlon  is  fundeaentally  different. 

This  rauld  noranlly  be  the  ense  for  the  degree  of  freedon;  tbnt  Is.  fomnrd 

speed,  not  distance  flora.  Is  the  ccntroll^  varieble.  Likeelse,  If  e  heeding 
refsreiice  oalj  is  uoSu  fot  uatiaatioa  (e.g.  ragnetic  (^nvaas).  then  u^  ,  not  x^  , 
heconen  the  rawkni  ootpnt.  In  sack  a  eaita  the  dlaplaceaant  In  n  givra  direction  is  , 
the  integral  of  the  correBpoodlng  randoa  velocity  co^onent.  l.e. 


If  m  consider  a  very  large  nuaber  of  flight  paths  througb  the  turbulent  field,  end 
take  an  ansaable  avanga.  denoted  hy  <  >  ,  than 


=  / 


<n|>dt  =  0 


ninee  the  anseable  average  la  equal  to  the  space  average.  The  aean  square  eo>ordlnate. 
hoeever,  does  not  vanish: 


t*(t)  =  Jj  Uj(a)Uj(/5)dad/9 


The  ansaubla  average  (average  over  aany  flights  at  givra  tlas  t)  is 


rj(t)>  =  <Uj(a)Uj(>8)>dad^ 


But  the  Man  product-  is  knoen  fre*  the  autocorrelation: 


R(a  -  /3)  = 


<Ui{a)mC8)> 
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Tlj«r«for« 


=1// 


R(a  -  ;8>dad/S 


liit«crml  can  be  sboen  to  hare  tbe  value 


R(T)<tr  -  2 


/ 


j 

» 


TR(T)dT 


<M) 


(67> 


At  larte  values  of  t  tlie  second  ten  becoaes  negligible,  and  ue  have  the  final 
result 

<ij>  =  2tt|At  (68) 

ulMru  A  =  /*8Cr)dr  is  the  area  under  tbe  antocorrelatlon  curve  of  Uj(t)  .  The 
latter  is  directly  related  to  tbe  ontpiit  poser  spectrua  of  u^  ,  and  can  be  ealeolated 
fra  It,  l.e. 


*<■»■)  -  J  («») 

:L  " 

ehere  6b^H|(c»)  la  the  poser  spectral  density  of  Uj  .  Hm  slgnlficence  of  tbe 
result  tlTSh  In  Bqwtlon  68  la  that  the  r.a.8.  disperaloi  (<x*>'i^  varies  ns  A  . 
this  is  tbe  sane  result  as  in  tbe  classical  problea  of  the  ’randan  salk* .  Ibus 
tbe  probable  error  in  tbe  lateral  position  of  a  oaspnas>CQatrollsd  fllidit  path 
increases  sltb  tbe  sgoare  root  of  tbe  tine,  or  distance  floua.  The  sene  scald  be 
true  of  tbe  distance  flos  Itself  in  a  speed*cQntrolled  flight.  Hosever.  since  an 
altitude  reference  is  alaost  invariably  used  la  tbe  fllibt  of  airplanes,  tbe  probable 
error  in  tbs  height  reneins  constant  sltb  tine.  The  dispersion  of  entirely  ungulded 
bodies,  e.g.  ballistic  aissiles,  would  vary  as  /t  in  all  three  co-ordinates. 


ACENt«Ln«KNT 

This  report  cougirises  the  srltten  version  of  an  Invited  paper  presoited  at  a 
neeting  of  tbe  A6MD  Flight  Hechanlcs  ftnel,  Brussels.  April  10-14.  1961.  Hiante 
are  due  to  Professor  E.D.  Poppleton  for  constructive  crltlcissi;  and  to  tbe  Nsticiinl 
Research  Council  of  Canada  for  sm^lylng  the  necessary  travel  funds. 
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and  Mil  Spec  8666  Requirements.  (Reproduced  from  Ref.  10.  Douglas 

Aircraft  Co. ) 
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(b)  Aa  distribution  associated  with  for  a  sweptback  wing 


Variations  of  angle  of  attack  along  span  of  a  high>aspect-ratio  wing 
associated  with  the, v-derivatives 


DISCUSSION 


O.E.  Michae  isen  ((^nada):  I  believe  Professor  Etkln  mentioned  that  the  difference 
between  the  simple  one-dimensional  gust  and  the  random  turbulence  may  be  important, 
particularly  with  regard  to  v/spit.  aircraft.  I  would  believe  that  this  indeed  is 
the  case  since  the  changes  in  inflow  angles  to  the  aircraft  due  to  gusts  are  of  the 
same  order  as  the  basic  flow  angle  la  the  steady-state  condition.  In  addition,  an  . 
amplification  of  the  turbulence  effects  can  occur  as  a  result  of  the  induced  changes 
in  the  slipstream  flow  angle  to  a  tilt,  or  deflected  slipstream  wing,  I  would 
appreciate  it  if  Professor  Etkin  would  comment  upon  this. 

Reply  by  Author:  I  agree  that  the  behaviour  of  vehicles  in  hovering  and  low-speed 
flight  subjected  to  low-level  turbulence  is  Important  to  understand  in  connection 
with  V/STOL  and  U/TUL  aircraft.  As  Mr.  Michaelsen  points  out,  the  Inflow  angles  may 
be  so  large  that  non-linear  aerodynamics  Is  involved.  This  will  certainly  present 
a  serious  difficulty.  It  may  also  be  expected  that  all  three  turbulence  components 
will  be  of  comparable  Importance,  and  that  the  correlations  associated  with  turbulent 
shear  flow  will  have  to  be  taken  into  consideration.  I  think  we  are  a 

long  way  yet  from  fully  understanding  the  behaviour  of  all  kinds  of  airborne  vehicles 
In  atmospheric  turbulence. 
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